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ABSTRACT 

Babcock-Leighton type solar dynamo models with single-celled meridional 
circnlation are snccessfnl in reprodncing many solar cycle featnres. Recent ob¬ 
servations and theoretical models of meridional circnlation do not indicate a 
single-celled flow pattern. We examine the role of complex mnlti-cellular circn¬ 
lation patterns in a Babcock-Leighton solar dynamo in advection- and diffnsion- 
dominated regimes. We show from simnlations that presence of a weak, second, 
high-latitnde reverse cell speeds np the cycle and slightly enhances the pole- 
ward branch in bntterfly diagram, whereas the presence of a second cell in depth 
reverses the tilt of bntterfly wing to an anti-solar type. A bntterfly diagram 
constrncted from middle of convection zone yields a solar-like pattern, bnt this 
may be difhcnlt to realize in the Snn because of magnetic buoyancy effects. Each 
of the above cases behaves similarly in higher and lower magnetic diffusivity 
regimes. However, our dynamo with a meridional circulation containing four cells 
in latitude behaves distinctly differently in the two regimes, producing solar-like 
butterfly diagrams with fast cycles in the higher diffusivity regime, and complex 
branches in butterfly diagrams in the lower diffusivity regime. We also hnd that 
dynamo solutions for a four-celled pattern, two in radius and two in latitude, 
prefer to quickly relax to quadrupolar parity if the bottom flow-speed is strong 
enough, of similar order of magnitude as the surface flow-speed. 

Subject headings: Sun: activity - Sun: magnetic helds - Sun: interior - Sun: 
helioseismology - Sun: photosphere 


2 


INTRODUCTION 


Over the past two decades Babcock-Leighton type OBabcocklllQGlI: lLeightonlll964119691) 
solar dynamo models operating with single celled meridional circnlation have been snccessfnl 
in reprodncing many solar cycle featnres, inclnding eqnatorward migration of s nnspot belts. 


1991; 

Ghoudhuri et a 

. 1995: 

DurnevI 

1995 

Diknati & Gharbonneau 

1999: 

Kiiker et al 

2991; 

Bonanno et al.l 

2992: 

Mandv & Ghoudhuri 

299 1|; 

Guerrero & Munoz 

2994; 

Jouve et al 

2998). 


It was possible to calibrate these models for the Snn and so they were applied to pre- 
di ction of solar cycle amplitnde . No w we know that sol a r cyc le 24 amplitnde forecast 
of Dikpati. de Toma fc Gilman ( 2006 ): Dikpati fc Gilman ( 2006 ). namely a 30% to 50% 
stronger cycle 24 peak than the peak of cycle 23, may not be validated. One of the reasons 
is that the assnmption of a steady, single-celled meridional circulation in each hemisphere 
may be oversimplihed for the Sun. Both observations and models indicate that there may 
be more than one cell in either depth or latitude, or both, in each hemisphere, at least at 
some times. 


Using time-distance helioseismology the most recent observatio ns from SDO/HMI data 
infer meridional circulation with two cells in depth (IZhao et al.ll2013[) . Rin g-diagram analysis 
from GONG data g i ves poleward surface flow up to about ~ 60° latitude (IHaber et al.ll2002 


Basu fc Antiall201Cll: iKomm et ah 1120131 ). whereas Doppler measurements from MWO data. 


which can measure this flow at higher latit udes, show a h igh-latitude, reverse flow associated 


with the primary poleward surface flow (lUlrich 1120101). Using a v ery long-term GONG 


database and applying time-distance technique iKholikov et al (120141) have found signatures 
of eqnatorward return-flow in the lower half of the conv ection zone, i ndica ting a long deep 
one cell flow-pattern. A p-mode perturbation analysis bv ISchad et al.l (120131) yields four cells 
in latitude, each going down to about 0.8R. Thus observations do not yet give us a unique 
answer about the Sun’s meridional circulation pattern. 


Models produce more complex flow patterns that vary from model to model. For ex¬ 
ample, mean-field models produce a long, counterclockwise primary cell often associated 
with a weak, reverse cell at high latitudes, both extending down to th e bottom of con- 


vection zone ( 

Rudiger 

1989; 

Kitchatinov & Riidiger 

2995; 

Diknati 

2914) 

, whereas full 3D 

convection simulations 

produce multiple cells in latitude and depth 

Guerrero et al. 

2913; 

Featherstone & Miesch 

2915). Recentlv 

Diknati. Anderson & Mitra 

(29U 

) have shown the 


scope of application of Ensemble Kalman filter data assimilation in a flux-transport dynamo 
model for reconstructing the time-variation of the Sun’s meridional circulation, and in the 
future the combination of theory, observation and data assimilation can be implemented to 
derive the spatio-temporal pattern of the Sun’s meridional circulation. However, given the 
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lack of knowledge about the uniqueness of meridional flow from observations and models at 
present, it is necessary to consider all plausible meridional circulation patterns for the Sun, 
and explore their effects on a Babcock-Leighton solar dynamo model. We specihcally seek 
the answers to the following questions: (i) can a Babcock-Leighton dynamo operating with a 
multi-cellular meridional circulation pattern produce observed solar cycle features? (ii) Can 
such a model be calibrated for the Sun in advection and diffusion-dominated regimes? 


In the present paper, our aim is to study the effects of various plausible multi-cellular 
meridional flow patterns on Babcock-Leighton solar dynamos operating in a full spherical 
shell of the convection zone. It is important to use inputs to the model, such as meridional 
circulation, that are as closely guided by available observations as possible, to get the best 
possible model dynamo properties that can be compared with properties of observed solar 
cycles. Guided by the observational and modeling results cited above, we choose hve circu¬ 
lation patterns. These include in each hemisphere (a) a single cell with a poleward flow at 
the surface; (b) a long primary cell from the equator to about 60° latitude, associated with 
a second, reversed cell at high latitudes; (c) two cells in depth; (d) two cells in depth and 
two in latitude; and (e) four cells in latitude. 


Flux transport dynamo simulations have been done for some of these cases, but before 
the latest helioseismic observations of meridional circulation was available, so the cho ices of 
circula tion patterns was less closely aligned to these observations. For example, [.Bonanno, Elstner, Belveder( 
020051) simulated a flux-transport dynamo operating with two cells in latitude having similar 
latitudinal extent and amplitude. The observations clearly show the low latitude cell always 
reaches to at least 60° and the second cell beyond it is quite weak compared to the primary 
cell (our case (b)). These differences lead to substantial differences in butterfly diagram, 
generally in the direction of poorer agreement with solar cycle observations. 


Jouve fc Brun 


fl2007l) have explored a Babcock-Leighton flux-transport dynamo with 
meridional circulation patterns equivalent to our cases (a), (c) and (e). However, there 
are a number of differences between our cases and theirs, which we explain below. Their 


motivation in large part was to try to improve on the success of solar flux transport dynamo 
models with single celled meridional circulation in simulating solar cycle features. However, 
they generally found that two- and four-celled circulations led to butterfly diagrams and 
other characteristics that were less like the Sun rather than closer to it. In this paper we use 
the single celled dynamo results as the reference case, and focus on which other circulation 
patterns suggested by most recent observations could also do as well as the reference case, 
and which give results that diverge signihcantly from the solar observations. 


In iJouve &: Brun (120071) all four cells have similar amplitudes and latitudinal dimen¬ 


sions, again signihcantly different from recent solar observations. In addition, the stream- 
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function for their meridional circulation is physically unrealistic for the Sun, because it 
is computed using a density prohle that varies like 1/r across the convection zone. This 
means that the density at the bottom is only 50% higher than at the outer boundary, so 
the flow is almost incompressible. The effect is to make the deep circulation cells much 
more nearly equal in amplitude to the surface velocities than is likely to be the c ase in the 
Sun. We will use the same density prohle as taken in iDikpati fc Charbonneaul (1199911 . in 
which density is proportional to {Rq/t — 7 )™). This density prohle gives a density diher- 
ence between top and bottom that is similar to a polytrope for an adiabatically stratihed 
so lar/stellar convection zone, a s well as to typical prohles used in helioseismology such as 
in IChristensen-Dalsgaard et al.l (1199611 . For values of 7 ~ 1 and m = 1.5, the density near 
the top becomes much more like the Sun for a bottom density of 0.2 gmcm“^. With this 
density p rohle, the bottom cell becomes much smaller in amplitude than is the top cell, as 
found bv IZhao et al.l (120131) from SDO/HMI data. These diherences in cell amplitude also 
strongly ahect the butterhy diagram and surface poloidal helds. 


The inference of two cells in depth by IZhao et al.l (l2013l) have created enough interest to 
explore how a hux-transport dynamo behaves with such a how pattern. Being poleward at 
the base of the convection zone, this how transports spot-producing tachocline toroidal helds 
poleward. Our case (c) will address this issue. Our case (d) is another form of two cells in 
depth, but the cells do not go all the way to the pole, instead reverse beyond 60° latitude. We 
study the role of such a four-celled pattern (two cells in depth and two in latitude), because 
Zhao et al.l (120131) did not conhrm whether the poleward surface how continues to the pole or 
stops at high latitudes. On the other hand, surface Doppler measurements in dicate that the 
surface how is poleward up to a certain high latitude (~ 60°); beyond that, lUlrich (120101) 
found a reverse, equatorward how during some epochs, although not all the time. 


We examine which of the hve circulation patterns mentioned in (a-e) can produce cyclic 
features similar to the Sun. Furthermore we perform a systematic parameter survey to 
compare dynamo model simulations for all hve circulations using the same model and the 
same dynamo physics, and judge which models that use these circulation patterns can be 
calibrated to the Sun in dihusion and advection dominated regimes. After describing the 
model in the next section, we present our results in §3 and conclude in §4. 
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2. MODEL 


2.1. Dynamo equations. 


We use the spherical polar coordinates (r, 0,0), and assume the axisymmetry. The 
magnetic held, as the sum of a toroidal component (B,^) and a poloidal component (Bp), 
can be written as: 

B(r, 0, t) = B 0 (r, 0, t)e^ + V x [A(r, 0, t)e^] (1) 

where the toroidal component of the magnetic held B^(r, 0 ,t)e 0 and the vector potential 
A(r, 0, t)e 0 . The both components can be generated by a how. The large-scale how held 
U(r, 0) be expressed as the sum of diherential rotation (h2(r, 0)) and the meridional circulation 
(u(r, 0 ) = Ur{r, 6)er + ue{r, 6)ee): 

U(r, 0) = u(r, 0)-I-rsin0r2(r, 0)e,ji (2) 


as toroidal and poloidal parts of the total how held. 

The evolution of the large-scale magnetic held B according to: 


dB 

dt 


= V X (U X B — 77 V X B), 


( 3 ) 


The toroidal component becomes 


dB^ 1 
dt r 


d d 

— {Tn,B^) + —{ueB^) 


= rsin0(Bp ■ V)!! — V ?7 x V x B^e^ 


(4) 

\ sm 0 / 

where 77 (r) is the magnetic dihusivity. 

(9A 1 / 1 \ 

^ += A + S(r,0;B<^) (5) 

c/t r sm 0 \ r"' sm 0 / 

where Bp = V x (Ae<^). 

There is considerable uncertainty about what is the most realistic prohle of turbu¬ 
lent magnetic dihusivity with radius. Direct measurements of magnetic dihusivity as a 
function of depth are not possible yet. The mixing-length theory gives us a rough esti¬ 
mate of the supergranular dihusivity in the super granulation layer near the surface with 
a range (77super = 10^^ — 10^^cm^s“^). The magnetic dihusivity in the convective enve¬ 
lope of the Sun is dominated by its turbulent contribution, but below the convection zone 
there is much less turbulence, the core is stabler, so the dihusivity should be determined 
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essentially from the molecular contri bution in the stably stratified deep radiatiye interior 
(fpikpati. Gilman fc MacGregor 2006 ). We assume that the turbulence governs the diffu- 
sivity in the convecti on zone and is signifi cantly reduced in the subadiabatically stratified 
radiative zone below flDikpati et ah 1120041) . The diffusivity profile can be written as: 


/ N , Vt 

h(r) = hcore + Y 


1 + erf 


r - rs 


+ 


^super 


1 + erf 


r - rg 


( 6 ) 


The parameters in the diffusivity profile: = 3 • 10^^ cm^s“^ is the turbulent diffusiv¬ 

ity, rjcove = 10® cm®s“^ is the core diffusivity, rjsupev = 3 • 10^® cm^s“^ is the supergranular 
diffusivity, rg = 0.7 R©, dg = 0.0125 Rq, rg = 0.9562 Rq, dg = 0.025 Rq. 


A two-step profile, which leayes out the supergranule layer, will have impact on the 
surface poloidal fields that are used to create the poloidal part of the butterfly diagram. 

Our flux-transport dynamo can be driven by both a tachocline a-effect and Babcock- 
Leighton type surface poloidal source. The decay of tilted bipolar active regions produce 
poloidal fields near the surface, but perhaps that cannot be the sole driver for the large- 
scale solar dynamo. A Babcock-Leighton dynamo is not a self-excited dynamo, and hence 
cannot come back to a normally cycling dynamo if it goes to a grand minimum state. Flux- 
transport solar dynamo models with both Babcock-Leighton type surface poloidal source ari d 
convection zone a-effect have been used (jPassos et al.ll20l4 iHazra. Passos fc Nandv 11201411 . 
in the context of reviving a Babcock-Leighton dynamo from a grand minima. Furthermore, 
it is well known that a flux transport dynamo driven solely by a Babcock-Leighton poloidal 
source cannot susta in its antisymme tric magnetic field about the equator, as inferred from 
Hale’s polarity rule flHale et al.lll919l) . but slowly drifts to symmetric mag netic field solutions. 


This i s true even in the case of a 3D Babcock-Leighton dynamo (see, e.g. iMiesch &: Dikpati 


(1201411 1. Inclusion of a tachocline instability-driven a-effect is demonstrated to be one of 
the solutions to stop such a drift to a nonsolar-like, quadrupolar parity in a pure Babcock- 
Leighton dynamo models flDikpati fc GilmanI 1200 ll: iBonanno et al.ll2002l) . In the present 
paper, since one of our goals is to investigate the parity produced by a Babcock-Leighton 
dynamo model operating with different multi-cellular meridional circulation patterns, we 
assumed zero tachocline a-effect in all cases. We use the following expressions respectively 
for Babcock-Leighton surface source: 


S(r, 9, B^) — {SBL(r 7 Stachocline(^5 0)} 


1 + 


B^(r,0,t) 


Br 


-1 


(7) 


*S'tachocline(^) 0 


( 8 ) 
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’BLl 


Si 




sin 9 cos 9 

2 

L V d4 


V 4 h 



X _l_ e7i (7r/4-e) 


( 9 ) 


for 6 < TT12 and for 6 > -n 12 
SBL(r,0) 


Si 




sin 9 cos 9 

2 

V d4 h 


V df. h 



X _|_ eT'l(^'-37r/4) 


( 10 ) 


The parameters used in ([9]) and ffTOj) are: Si = 2.5 ms = 0.95 R©, rs = 0.9875 R©, 
d4 = dg = 0.0125 Rq, 7i = 30.0, (3 = 70.0, and Bq = 10 kG. 

The solar internal rotation prohle includes primarily latitudinal shear in the convection 
zone, as found by two-dimensional helios e ismic inversions. A solar-like internal differential 
rotation prohle fjpikpati fc Charbonnean (1999)) is given by 


f2(r, 6^) — Gc + 2 


1 -|- erf ( 2 


r — rr 


di 




(11) 


where 

a = + a2Cos^ 9 + a4Cos'^ 9 (12) 

is the surface latitudinal differential r otation. Values wer e chose n to closely resemble the 
best ht to the helioseismic solution of ICharbonneau et al.l fjl998[) . The angular velocity of 
the rigidly rotating core is VLc/2ti = 432.8nHz. GEq/2vr = 460.7nHz is the rotation rate at 
the equator. The other parameters are set to be a2/27r = 62.69 nHz, and a4/27r = 67.13 nHz. 
The rc = 0.7 R© indicates the central radius of the tachocline thickness di = 0.025 R©. 

The differential rotation prohle prescribed by three terms, as shown in Equation (12), is 
well-formulated up to 60°, but to adequately ht the measured rotation rate including more 
poleward latitudes requires additional terms flSchou et al.lll998l: iDikpati et ah II2002I) . We 
will discuss how the implementation of such diherential rotation prohle ahects our results. 


2.2. Stream function 


The circulation is represented in the spherical shell by the stream function of iForgacs-Daika fc Petrovav 


(12OO2I) for both hemispheres. To study the ehects of the meridional circulation, we use a 
simple and easily adjustable spatial structure, with realistic amplitude. The components of 
the meridional circulation can be written as 

1 d 


Ur{r, 9) = 


p{r) ■ r'^ sm.9 09 




(13) 













































ue{r,0) = 


d 


^{r,e) 


p{r)r sin 6 dr 

Note that the spherical polar geometrical factor (r sin0) is absorbed in ^{r, 6 ). 


(14) 


(o) (b) (c) (d) (e) 




Fig. 1.— Streamlines of the meridional circulation as prescribed in flTHl) . 


We reproduce the stream function of iForgacs-Daika fc Petrovavi fj2002l) . which has the 
following form: 

^'(r, 6 ) = sin^ ni 6 cos n 26 (15) 

where the form of the given function ■^(r) specihes the flow. 


V’(r) = V’o sin 


kvr (r — r,, 




O ■‘■me 


exp 


(r - ro)" 

F2 


(16) 


where V’o sets the amplitude of the meridional circulation. Two parameters play the pri¬ 
mary role in setting the meridional circulation for the simulations. We £x the number of 
cells in latitude by changing the parameter n and the number of cells in radius by chang¬ 
ing the parameter k. The rg and F are geometric parameters, rg = (R© — rjnc)/30cm and 
F = 6.4 ■ 10^° cm. rmc is the radius to which the meridional flow penetrates from the base of 
the convective zone. Given the observed depth of the tachocline, we set rmn = 4.79 ■ lO^*’ cm. 
We use the same density prohle (p(r)) as used bv lDikpati &: Charbonneaul fjl999l) . 


With all parameters specihed, we can compute the stream functions for the hve circu¬ 
lation patterns we use in the dynamo. Figure 1 displays the circulation patterns we will 
use for our analysis of the differences in dynamo behavior that arise from different circula¬ 
tion patterns. These include the reference case (a) a single cell in each hemisphere; (b) a 




































9 


long, primary cell associated with a weak reverse cell in high latitudes; (c) two cells with 
depth; (d) two cells in latitude and two in depth; and (e), four cells in latitude. There 
is some observational support for each of these patterns, as reviewed in the introduction. 
The parameters chosen in the hve cases are given respectively as: (a) /c = 1, ni = n2 = 1, 
To = —6.1771 X 10^^; (b) A: = 1, ni = 1, ^2 = 3, Tq = 3.519 x 10^^; (c) k = 2, ni = 1, n 2 = 1, 
To = 3.0937 X lO^i; (d) A: = 2, ni = 1, n2 = 3, To = -1.7625 x lO^i; (e) A: = 1, ni = 4, 
n2 = 4, and Tq = Psio^ + + ^Oc + in which we assign values for Tq^, Tq;,, Tq^, Tq^ 

for a certain range in 6 and zero elsewhere, namely Toa = 1.7043 x 10^^ for 0 < 0 < 22.5°, 
Tofe = 4.1340 X 10^1 for 22.5° < 9 < 45°, Tq^ = -6.8972 x lO^i for 45° < 0 < 67.5°, 
Tod = —7.8002 X 10^^ for 67.5° < 6 < 90°. Unit of Tq is in c.g.s, i.e., cm^s“^ unit of p~^. 


2.3. Boundary conditions and solution method 


Equations (jl]) and (E]) are solved in a full-spherical shell (0 > 6 * > tt), extending radially 
from below the bottom of convection zone (r/RQ = 0.6) to the surface (r/RQ = 1). We use a 
4th order Runge-Kutta method with the central hnite difference scheme for r and 6 deriva¬ 


tives to solve the equations. Validity test of the scheme was performed by reproducing 


suits the benchmark cases of a Babcock-Leighton flux-t ransport dynamo flJouve et al 


he re- 


20081 ). 


We use the boundary conditions as used previously by Dikpati fc Charbonneau ~ fjl999 '). We 


set = 0 and A = 0 at the bottom boundary {r/R q = 0.6). Note that the bottom 
boundary condition for the toroidal field i s modi fied from the condition for a perfectly con¬ 


ducting core (see e.g. iDikpati fc GilmanI (120011) for physical explanation). Both B,^ = 0 


and A = 0 are also set to zero along the polar axis (6* = 0 or vr) to ensure physicality. At 
the surface, the toroidal fields are set to zero, and we demand that the poloidal held-lines 
match smoothly with the p otential field solutions of the free space above the surface (see 
Dikpati fc Choudhuril (119941) ). The vector potential A satishes the following equation above 
the surface: 

V" - \A = 0. (17) 


(18) 


sin^ 9 ^ 

A general solution for this equation can be written as: 

A{r > Rq, 9,t) = Y^ ^^P,((cos 9) 




(2n + 1)R^ 


n+l 

© 


t/2 


n{n + 1) 


A{r = RQ,9,t)Pl^{cos9) sm9d9 


(19) 
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where P^(cos 6 ') is the associated Legendre polynomial. The derivative of A at the solar 
snrface 

ir=RQ) n 

3. RESULTS 
3.1. Time-latitude diagrams for multi-cell flow 

First we establish a reference case dynamo solntion, to which other individual solutions 
can be compared, to see what changes are created by changing the meridional circulation 
pattern. We choose for reference the frequently used single cell meridional circulation that 
has poleward flow near the outer boundary, and return flow at the base of the convection 
zone (see Figure 1(a) and Figure 2(a)). This flow penetrates slightly below the base of the 
convection zone, given by r = 0.7Rq. The thin black dashed semicircular arc represents 
the location of the center of the tachocline. Two frames, (b) and (c) on the right panel 
of Figure [21 show the time-latitude diagrams respectively of toroidal held at the bottom of 
convection zone (B(^|i.=o. 7 R 0 ) and surface radial held (Br|r=RQ). To get these results we took 
a maximum surface flow speed of 15ms“^, poloidal source amplitude of Si = 3.0 ms“^ and 
turbulent diffusivity = 3 ■ 10 ^^ cm^s“^. 

Figure [3] shows how the poloidal and toroidal fields evolve through a sunspot cycle for 
the entire meridional cross-section. When the toroidal field is strong, it is confined to the 
lower layers of the domain, and its peak is clearly migrating toward the equator, along with 
the poloidal field lines that are sheared by differential rotation to produce it. We expect 
sunspots to emerge from at or near the latitude of maximum in the toroidal field. The 
reversal of polar helds is seen here in the interval between 5 and 6.25 years, at which time 
the toroidal is a maximum with its peak near 20 °. 

Babcock-Leighton flux transport dynamo models with one-celled meridional circulation 
in each hemisphere can reproduce many features of the solar cycle. These include (i) the 
equatorward migration of toroidal flux at lower latitudes; (ii) the 11 -year sunspot cycle; 
(hi) the observed phase-shift between poloidal and toroidal components; (iv) the short rise 
of toroidal field to maximum followed by the long decline to minimum - in the reference 
case, the ascending phase is 16.77% of the whole cycle; (v) peak tachocline toroidal fields are 



dA 

dr 
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0 20 40 60 80 

time (year) 


Fig. 2.— Panel (a) displays the streamlines for single-celled meridional circnlation in 
each hemisphere; bine-violet represents connterclockwise flow, green-orange clockwise flow. 
Panel(b) shows the time evolntion of the tachocline toroidal held; panel (c) the same for 
snrface radial helds. For tachocline toroidal helds, color-hlled contonr levels are 3kG; the 
highest tachocline toroidal held is ~ 37kG (yellow/violet). The maximnm valne of the radial 
helds is ±204 G, occnrring near the poles. 


Belucz &: Dikpatil 1201311 . This solntion also reprodnces the observed phase shift between 


the surface poloidal held and the toroidal held at the tachocline; the poloidal held polarity 
changes from positive to negative when the toroidal held is near maximum and its polarity 
is negative. 


In the next simulations we study how the characteristic features of butterhy-diagram 
change when the meridional circulation cell contains a second, high-latitude, reverse cell. 
The form of this meridional circulation is shown in Figure 1(b) and Figure 0^. The peak 
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Fig. 3.— Evolution of toroidal and poloidal fields with one-cell meridional circulation (Figure 
[T](a)). The hlled contours show toroidal helds, yellow/red colors denoting positive (into the 
plane of the paper) and blue/purple negative (out of the plane). Red and blue contours 
respectively denote positive (clockwise) and negative (anticlockwise) poloidal heldlines. 

flow-speed of the primary cell is still 15ms“^, poleward at the surface, but the latitude of 
this peak is slightly lower, at 25.3°. The peak flow speed of the secondary cell is 3ms“^, 
equatorward at the surface. The boundary between cells is near 61° latitude. 

The right panel of Figure 0] shows the time-latitude diagrams of B<^|r=o. 7 R 0 in panel 
(b) and Bi.|r=R 0 in panel (c). Not surprisingly, due to the effect of the second cell, a more 
pronounced poleward branch can be seen in the butterfly diagram of toroidal held in panel 
(b) compared to that in Figure [2l The sunspot cycle length (i.e. half magnetic cycle) is just 
8.3 years, due to the shorter primary conveyor belt. The strength of toroidal held is similar 
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0 20 40 60 80 

time (year) 


Fig. 4.— Same as in Figure |2] but the meridional flow pattern has a high-latitude second 
reverse flow-cell (Figure [T]d). The highest tachocline toroidal held is similar to that in Figure 
[21 ~ 37kG (yellow/violet). The maximum value of the radial fields is ±89 G, at about 50° 


to that of reference case, 37 kG. Gomparing the time-latitude diagrams of Figure [21(c) and 
[11(c), we And that the polar held peaks around 50° latitude instead of peaking near the pole 
as in the case of a single cell in each hemisphere. This is due to flow convergence at 61° 
latitude instead of at the pole. The second cell also causes a delay in the polarity change by 
advecting polar fields away from the pole. 

The rise of the cycle from minimum to maximum in this case is slightly longer compared 
to that in the single cell case. This is probably because some of the poloidal flux advected 
to the bottom in between the primary and secondary cells is advected toward the poles, 
retarding the early production of the equatorward migrating branch of toroidal held there. 
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For the next simulation, we add a second, reversed meridional cell below the primary 
cell (see Figure [It). The two cells are similar in amplitude and radial extent. The results 
for this meridional flow are shown in Figure O 
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(b) Toroidal field at 0.7R (c) Radial field at LOR 



40 60 

time (year) 


Fig. 5.— Same as in Figure |2] but the meridional flow pattern has a second, reversed flow¬ 
cell below the primary cell (Figured}:;). The maximum tachocline toroidal held strength is 
~ 23kG (orange/violet). Maximum value of the radial fields is ±83G, near the poles. 


This dynamo solution is very different from the reference case. At all latitudes, the 
tachocline toroidal held is migrating with time toward the poles rather than the equator. 
The rate of this migration increases with latitude. This is due to the poleward flow at 
the bottom of the second cell. Thus, based on tachocline toroidal fields, this meridional 
circulation pattern produces an ’antisolar’ butterfly diagram. The toroidal held at the bottom 
of the domain is weaker than in the reference case, because polar fields are not brought down 
from the top, but instead are advected from lower latitudes near the bottom. 

If we plot instead the toroidal held near the middle depth, for example, 0.83125 R, 
shown in Figure 6, we get a more solar-like butterfly, with both poleward and equatorward 
branches. This is because at these depths both circulation cells have equatorward flow, so 


























- 16 - 



time (year) 


Fig. 6.— Same as in Figure [5] but the toroidal fields are shown from 0.83125 R. A solar-like 
butterfly diagram is obtained to equatorward flow there. 

they advect toroidal field toward the equator in lower latitudes. The relatively high speed 
total flow there also makes the dynamo period slightly shorter. In addition, since in this 
case the poleward flow near the outer boundary reaches to a shallower depth, less poloidal 
flux is advected toward to the pole. These weaker polar fields lead to weaker toroidal fields 
at all latitudes. 

The next dynamo simulation is for the case of two meridional cells in both latitude and 
depth in each hemisphere. The amplitudes of upper and lower cells is about the same. These 
results are shown in Figure [71 Panel (b) again shows the toroidal fields near the bottom of 
the convective zone. Panel (c) depicts the surface radial field. As in the case of two cells in 
radius, we again get an ’anti-solar’ butterfly diagram when we plot the toroidal field at the 
bottom. Again this is because the flow is toward the pole at the bottom of the stack of cells 
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time (year) 


Fig. 7.— Same as in Figure |2] but the meridional flow pattern has four flow-cells (Figure 
[T]i). The highest tachocline toroidal held is ~ 33kG (yellow/violet). The maximum value 
of the radial fields is ±100 G, at around 50°. 


at all latitudes except the highest. If we plot the toroidal held contours at mid-depth (not 
shown), we will again get a solar type butterfly, as seen in Figure 6. 

In the last simulation, we also have four cells, but these cells are located side by side, 
as seen in the Figure [Tfe. In this case, the fields are confined to lower latitudes because the 
multiple cells in latitude prevent poloidal held transport all the way to the poles as in the 
reference case. We get a solar-like time-latitude diagram up to about 22°; the toroidal held 
migration is equatorward. We do not see dynamo activity beyond about 50°. The cycle 
length is very short, 3.125 years, due to the very short conveyor belts represented by the two 
circulation cells closest to the equator. The strength of the toroidal held is just half of the 
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Fig. 8.— Same as in Figure |2] but the meridional flow pattern has four flow-cells (Figure 
d^). The highest tachocline toroidal held is ~ 15kG (red/violet). The maximum value of 
the radial fields is ±33 G, which occurs near 25°. 


reference case, because the dynamo is confined to the lower latitudes where the differential 
rotation is smaller so the production of toroidal held from a given poloidal held is smaller. 
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3.2. Time-latitude diagrams for multi-cell flow with lower maguetic diffusivity 



'b) Toroidal Held at Q.7R (c) Radial at LOR 



40 60 

time (year) 

Fig. 9.— Same as in Figure [2] but the turbulent diffusivity is 7 • 10^°cm^s“^. The highest 
tachocline toroidal field is 52.4kG (red/violet). The maximum value of the radial fields is 
±207 G near the poles. 


How dependent are the results described above on the particular choice of turbulent 
magnetic diffusivity? We address that question by repeating the dynamo simulations for the 
five meridional circulations using a lower value, 7 • 10^°cm^s“^. We selected this value of 
diffusivity {t]t = 7 • 10^°cm^s“^) to represent the low-diffusivity cases, and in the previous 
subsection §3.1 a diffusivity value of 3 ■ 10^^ cm^s“^ to denote the representative value for 
high diffusivities. However, these values are not unique; neighboring values can be considered 
also. 


The results for low diffusivity are displayed in Figures 9-13. Figure |9] displays the 
reference case with single circulation cell for the lower magnetic diffusivity. As in earlier 
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figures, panel (b) gives the toroidal field amplitude near the bottom, panel (c) the surface 
poloidal held. We use the same peak surface how speed (15ms“^) and poloidal source 
amplitude (si = 3.0 ms“^) as for the solution seen in Figure 2. 

In this lower dihusivity case, the model still reproduces the observed phase shift be¬ 
tween the surface poloidal held and the tachocline toroidal held. Due to the lower turbulent 
dihusivity, much stronger helds are produced. The highest tachocline toroidal held is 52 kG, 
and the maximum value of the radial helds is ±207 G near the poles. The period is also 
somewhat longer than in Figure 2, namely 12.7 years, as measured by the time between 
adjacent peaks in toroidal and poloidal held at the same latitude. But there is also sig- 
nihcant overlap between adjacent cycles, so the time between the high latitude peak of a 
new sign of toroidal held to its disappearance at low latitudes is more than 20 years. How 
much overlap there is in the Sun itself is unclear, since from observations we know where 
the tachocline toroidal held is only from the latitude of sunspots seen. The slight poleward 
migration of tachocline toroidal helds seen at high latitudes comes from the strong negative 
radial shear there overcoming the relatively weak equatorward advection of toroidal held by 
the meridional how there. 

Figure [10] shows the time-latitude diagrams of toroidal and radial helds in the two 
cell case fFigure [TOhi with 7 ■ 10^°cm^s“^ turbulent dihusivity. For this case we see that the 
high and low latitude branches of the butterhy diagram for tachocline toroidal held are about 
the same, despite the unequal amplitudes of the two circulation cells. Both the poleward 
meridional circulation at the bottom and the radial gradient of rotation in high latitudes at 
tachocline depths are contributing to this pattern, which is more pronounced than seen in 
Figure m for which the bottom meridional circulation is toward the equator at high latitudes. 
The sunspot cycle is shorter, about 10.7 years, due to the shorter conveyor belt that is the 
primary circulation cell. Here too, the induced toroidal field are stronger than in the high 
diffusivity solution with the same meridional circulation. The maximum value of the toroidal 
held is 78 kG; maximum radial held is 77 G. In this case we get substantially less overlap of 
adjacent cycles than in the low dihusivity reference case (Figure 9). 
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(b) Toroidal field at 0.7R (c) Rudial field at l.OR 
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Fig. 10.— Same as in Figure H] but the turbulent diffusivity is 7 ■ 10^°cm^s“^. The highest 
tachocline toroidal held is ~ 78. kG (red/violet). The maximum value of the radial helds is 
±77 G, occurring near 50°. 
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In the next simulation with low diffusivity, we turn on a second meridional cell below 
the primary cell, as for the case previously displayed in Figure |5l The results are shown in 
Figure [TT] The most prominent feature of the butterfly diagram is the longer sunspot cycle, 
about 50 years. What little migration of toroidal held there is, is toward the poles, leading 
to a slightly antisolar butterhy. Here again, as in the higher diffusivity case shown in Figure 
El toroidal held contours (not shown) at mid-depth of the convection zone, yields a butterhy 
diagram with both poleward and equatorward branches. 



0 50 100 150 

time {year) 


Fig. 11.— Same as in Figure El but the turbulent dihusivity is 7 • 10^°cm^s“^. The highest 
tachocline toroidal held is ~ 114kG (red/violet). The maximal value of the radial helds is 
±288 G, near the poles. 
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Fig. 12.— Same as in Figure [7] but the turbulent diffusivity is 7 ■ 10^°cm^s“^. The highest 
tachocline toroidal held is ~ 55 kG (red/violet). The maximum radial held is ±83 G, from 
about 50° to near the poles. 


The results of the next simulation is shown in Figure [I2i Here, with two strong cells 
in low and mid-latitudes, and two reversed, relatively weak cells in polar latitudes, we again 
get a predominantly antisolar butterhy but with a small domain of equatorward migration 
of surface poloidal helds in low latitudes. This must be caused by the equatorward how at 
mid-depth coupled with upward magnetic dihusion to the surface. Finally, here again we 
get a solar like butterhy for low and middle latitudes when we plot toroidal held contours at 
mid-depth in the dynamo domain. 
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Fig. 13.— Same as in Figure |8] but the turbulent diffusivity is 7 ■ lO^^cm^s The highest 
tachocline toroidal held is ~ 65 kG (red/violet). The maximum value of the radial helds is 
±27 G. 

In the last simulation, we show results for the four cell case for low diffusivity in Figure 
fT3l Due to the short conveyor belts, decreasing the turbulent diffusivity does not signihcantly 
change the cycle length. But the lower turbulent diffusivity has other effects. First, the helds 
are stronger, as we should expect. Second, we can see dynamo activity at the higher latitudes 
than in the high dihusivity case, though it is still low compared to that in low latitudes. The 
orientation of the wings of the time-latitude diagram at the diherent latitudes is determined 
by the latitudinal direction of the how near the bottom; equatorward how leads locally to 
migration toward the equator, and poleward how to migration toward the poles, as seen in 
Figure fT3b . 

From all the time-latitude plots shown above, we can see that the most solar-like di- 
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agrams are produced if there is a single primary circulation cell in each hemisphere, with 
possibly a weaker secondary, reversed cell in polar latitudes. Circulation with two cells in 
depth, or two cells in both latitude and depth, give solar-like butterflies only from toroidal 
helds at mid-depth, not the bottom. For these to be correct for the Sun, the toroidal helds 
at the bottom must not come to the surface because of their magnetic buoyancy or for any 
other reason, and a mechanism must exist that keeps mid-depth toroidal helds from ris¬ 
ing buoyantly too fast to be amplihed to produce sunspots. Neither requirement is easily 
satished using known MHD processes. 


3.3. Parameter survey 


We have shown how solutions from a solar hux transport dynamo model differ for 
different forms of meridional circulation. The solutions we obtained are all for the same 
differential rotation, since from helioseismic measurements that is relatively well known for 
the Sun. These were found for hxed meridional circulation amplitude, but that amplitude is 
not very constrained from observations, and it is bound to have time variations. In addition, 
there are uncertainties in the amplitude and form of the surface poloidal hux source as well 
as the turbulent magnetic dihusivity. It is of interest to know how basic characteristics of 
a simulated sunspot cycle diher for diherent values of these uncertain parameters. Three 
prominent features of simulated cycles to focus on are its period, amplitude, and shape - 
the times spent in ascending and descending phases. We focus here on the hrst two of these. 


We choose the range of meridional circulation amplitudes to cover the range observed, 
and somewhat beyond. The choice of surface poloidal source amplitu de range is guided 
by est imates of relation to photospheric magnetograms as discussed in iDikpati fc Gilman 
(120061) . While the choice of ranges of values for the meridional circulation amplitude and 
surface poloidal source is guided by observations, and hence adequately map the plausible 
parameter space for the Sun, the range of the turbulent diffusivity amplitudes is selected 
based on the values that produce a sustaining dynamo. Given the turbulent diffusivity is one 
of the biggest unknown ingredients in the solar interior, we performed extensive numerical 
experiments to choose the range of rjT so that the dynamo does not die due to too much 
diffusive decay, or does not produce unusually large magnetic held amplitudes due to too 
low diffusivity. 
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Fig. 14.— Dependence of simnlated snnspot cycle length on (a) amplitnde of meridional 
circulation, (b) turbulent diffusivity and (c) amplitude of poloidal source term for the five 
circulation patterns used. The case of two cells in both latitude and depth is shown only 
for circulation amplitudes > 15ms“^. Below that amplitude, the solutions are quadrupolar 
rather than dipolar. 
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Figure [H] displays the variation of cycle period with circulation amplitude (panel (a)), 
turbulent diffusivity (panel (b)) and poloidal source amplitude (panel (c)). We see from panel 
(a) that as the circulation amplitude is increased, in almost all cases the period declines. 
This is to be expected, because in all cases, unless diffusion dominates, the period is set 
by the speed of the conveyor belt. The primary exception we see is, that for low speeds, 
decreasing the circulation in the single cell case (and to much lesser degree, the case with 
a second weak cell at high latitudes) leads to a decrease in period. This happens because, 
while the circulation is decreased, turbulent diffusion starts to ’short circuit’ the conveyor 
belt, since some poloidal flux is diffused toward the bottom from the top before it reaches 
polar latitudes. 

This short circuiting effect is even more evident in panel (b), where we have plotted 
cycle length versus turbulent diffusivity. For the same circulation amplitude, the solutions 
become more diffusion affected to the right in the figure. The periods decline, in some cases 
by factors of five or more. This result shows that to have a flux transport dynamo calibrated 
to the observed sunspot cycle period requires careful choice of the turbulent diffusivity, no 
matter what circulation pattern is assumed. By contrast, panel (c) shows that the cycle 
period is almost independent of the amplitude of the surface poloidal source. This is also 
expected, because the dynamo is nearly linear. Changing the poloidal source amplitude 
should change primarily the peak amplitude of the cycle, as we shall see below. 

Figure [T5] shows how the maximum toroidal held varies with meridional circulation 
(panel (a), turbulent diffusivity (panel (b)) and amplitude of the poloidal source term (panel 
(c)) for the hve circulation patterns. From panel (c) we see that, as we should expect, rais¬ 
ing the amplitude of the poloidal source raises the peak toroidal held amplitude. Because of 
the nonlinear quenching of the source term internal to the model, however, the amplitudes 
are beginning to approach asymptotic limits. There are also signihcant diherences in the 
efficiency of diherent circulation patterns; two cells in depth and four cells in latitude both 
generate much less toroidal held than the other patterns, which give almost the same ampli¬ 
tudes. In the four cell case, this is because the shearing of poloidal held is largely conhned 
to low latitudes, where the latitudinal rotation gradient is weakest, by the short latitudinal 
extent of the conveyor belt, so less toroidal held is generated. In the case of two cells in 
depth, less of the poloidal hux gets to the bottom where the radial shear is strongest, again 
reducing the dynamo’s ability to amplify toroidal held. 
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Fig. 15.— Dependence of maximum toroidal field on (a) meridional circulation amplitude, 
(b) turbulent diffusivity and (c) amplitude of poloidal source term, for the five circulation 
patterns used. The sudden drop in cycle amplitude for the case of two cells in latitude and 
depth is due to the solution switching from dipolar to quadrupolar at a peak circulation 
amplitude of 15ms“^. 
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In panel (b) the result that peak toroidal held declines with increasing magnetic dif- 
fusivity is due to the fact that there is more dissipation in the system that the induction 
processes in the dynamo have to overcome. Here too we see that, for all diffusivities, conhg- 
urations with four cells in latitude and two cells in depth generate less toroidal held for the 
same assumed dihusivity. 

From panel (a) we see that the variation in peak toroidal held with changes in circulation 
amplitude are more complex. With two cells in depth the peak amplitude changes relatively 
little with circulation amplitude. With two cells in both latitude and depth, we see a similar 
result, except near circulation amplitude 15ms“^ where the amplitude drops by more than 
50% when the magnetic held conhguration switches from dipolar to quadrupolar. In both 
cases, with the relatively hne scale meridional circulation pattern, the solutions are in the 
dihusivity dominated regime, so the amplitude does not change. Meridional circulation acts 
mainly as a transporter of toroidal and poloidal hux, rather than as an ampliher. The 
faster the toroidal held is transported, the weaker is the resulting held amplitude, because 
the toroidal held peak is displaced before it amplihed as much as it could be if it were 
stationary. 


The other three cases each show a largest value of toroidal held amplitude at a meridional 
circulation amplitude between about 8ms“^ and 15ms“^. At these speeds there is an optimum 
balance between amplihcation of toroidal held by diherential rotation shearing, dihusive 
decay, and meridional transport of toroidal and poloidal hux. In each case, for larger than 
optimum circulation, the toroidal helds are moved in latitude and/or depth too fast to be as 
fully amplihed as it would have been if moving more slowly; at less than optimum advection 
rate, more time is allowed at a given latitude and depth for the toroidal and poloidal helds 
to decay due to dihusion. 


Jouve fc Rrnn (120071 1 computed a power law relationship among the various parameters 


varied for the case of two cells in latitude and two in radius. However, for variation in how 
structures in terms of not only number of cells in latitude and depth, but also the node 
locations in latitude and depth, it will not be possible to reach a unique answer for the power 
law relationship. Thus we presented here the curves as function of how-speed, dihusivity 
and poloidal source amplitudes. Furthermore, some of these how patterns do not produce 
solar-like features, and hence power law relationships may not be meaningful. 
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3.4. Parity issue 


In all the cases we have presented so far, we have found dipolar parity during the 
simulation time of interest, i.e. about up to 500 years. We know that the growth rates 
of quadrupolar parity solutions are slightly higher than the dipolar parity in a Babcock- 


simulations run for 

more 

than 2000 years (! 

Dikoati & Gilman 
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Bonanno et ah 

2002; 

Hotta & Yokovama 

2010; 
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2014 

). However, 

Jouve & Brun 

(2007 

have 


shown that this switching from dipolar to quadrupolar parity is very fast in the case of a 
four-celled meridional circulation that consists of two cells in latitude and two in depth. In 
our four-celled case consisting of two cells in latitude and two in depth (see Figure [7]), the 
parity change did not occur so quickly. This is because the ratio of polewa rd surface flow- 
speed of the top cell to that of the bottom cell was too high, ~ 50, whereas in ijouve fc Brun 
(120071) that ratio was ~ 6. In order to investigate the fast change of parity in the four-celled 
case, we consider a ratio of poleward surface flow-speed to poleward bottom flow-speed to 
be ~ 5 and simulate that case and present our results in Figures 16-19. 


Even though the four-celled case with two cells in latitude and two in depth produce 
least solar-like solutions, it is worth examining the parity issue, because such a four-celled 
profile could appear intermittently in one or both hemispheres of the Sun. Fast change in 
parity in one hemisphere with respect to the other could lead to a large phase shift between 
North and South cycles, as observed in recent cycle. 


In Figure 16, as in earlier hgures panel (b) shows the toroidal helds near the bottom of 
convective zone. Panel (c) depicts the surface radial held. We see here a radical change from 
the earlier examples. The toroidal and poloidal held patterns are now symmetric rather than 
antisymmetric about the equator. In other words, we have found quadrupolar type rather 
than dipole type parity. This diherence in parity about the equator develops in just a few 
cycles, so the system in this case has a strong preference for quadrupole parity. 


Jouve fc Brun (1200711 demonstrated, in their simulations of dynamo with four-celled 


(two in latitude and two in depth) meridional circulation, for what combinations of merid¬ 
ional how speed and turbulent magnetic dihusivity quadrupolar structur e is favored. In the 
case o f a meridional circulation consisting of single cell in each hemisphere, iDikpati &: Gilman 
(120011) reasoned that slow switching to quadrupolar from dipolar parity occurs when bottom 
poloidal helds become weak enough after a long traversal via the conveyor-belt from surface 
to the bottom, and hence cannot connect with their opposite-hemisphere counterparts about 
the equator. A global statement would be that the quadrupolar mode is selected when, for 
a particular meridional circulation, it is dissipated at a substantially lower rate than is the 
dipolar mode. In other words, the growth rate for a quadrupolar mode is higher than for a 
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dipolar mode. Since a circulation pattern that has two cells each in both latitude and depth 
is inherently more complicated than a simple, single-celled pattern, the poloidal and toroidal 
flux of opposite signs and different amplitude are converged together in more places away 
from the equator. So the dynamo favors quadrupolar symmetry, unless the upper and lower 
cells are very unequal in amplitude. 
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Fig. 16.— Same as in Figure [7] but the bottom cell of this four-celled meridional circulation 
pattern has about ten times stronger poleward flow compared to that in [7^. The highest 
toroidal held amplitude is about ~ 27kG (yellow/violet). The maximum value of the radial 
fields is ±150 G, near 50° 


We And also in the low diffusivity case that having two cells in both latitude and radius 
leads to a quadrupolar solution, if the ratio of surface-flow to bottom-flow is 5:1 instead 
of 50:1 (we do not produce here the time-latitude diagram for the low diffusivity case). 
Therefore this fast switching from dipolar to quadrupolar parity is a typical phenomena in 
























the four-celled meridional circulation with a much stronger bottom-flow of similar order of 
magnitude as the surface-flow, irrespective of turbulent diffusivity value. 

Figures 17 and 18 give details of the actual transition from dipolar to quadrupolar 
symmetry. In Figure 17, panels (a) and (c) we see by eye that the transition appears to 
occur in about 120 years (16-18 sunspot cycles) in the tachocline toroidal held and surface 
poloidal held respectively. Frames (b) and (d) show hner detail for each in the middle 
of this transition. The simulation was started from a previous dipole simulation, which 
inevitably has some slight departures from dipole symmetry at the truncation error level. 
The quadrupolar symmetry is so strongly preferred for the parameters chosen that even these 
small diherences are enough to start the process of symmetry switching. 

What we see in Figure 17 is a very simple process in which switching occurs by the 
northern hemisphere developing a phase lag relative to the southern hemisphere, which grows 
until the North lags by a sunspot or half magnetic cycle, with very little change in pattern 
in each hemisphere. What is actually happening is that the periods of both hemispheres 
are getting longer, from about 6.Syr to 8.Syr, but the North reaches the longer period faster 
than the South, so the South gains on it in phase. Presumably in another simulation with 
different truncation errors, it could be the northern hemisphere that lags, ending up with 
the same hnal state. 

In Figure 18, we show meridional cross-sections of both toroidal held (blue/yellow shad¬ 
ing) and poloidal held (solid black and dashed red lines) in 1-year intervals in the middle of 
the transition. We can see particularly that the peak toroidal helds are moving up in the 
domain to a mid-depth where the how in low latitudes is toward the equator. In addition, 
the polar held reversals are evolving to a state in which they go from positive to negative at 
both poles at the same time. 
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Fig. 17.— Panels (a) and (c) show time-latitnde diagrams for tachocline toroidal fields and 
snrface radial fields respectively; evolntion of parity from dipolar to qnadrnpolar is shown in 
enlarged form in panels (b) and (d). 
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Fig. 18.— Twelve snapshots of toroidal fields in orange/violet color-filled contours and 
poloidal fields in red (positive) and blue (negative) contours show the evolution of helds 
during parity change. 
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We can quantify how long the transition takes by measuring the difference in cycle 
periods between South (P5) and North (P/v) as a function of time. The normalized difference 
in cycle periods between South and North (P^ — Pj^)/[{Ps + P/v)/2] is shown as function of 
time in Figure 19. We see that during this transition, the North adjusts to the longer period 
faster than the South. The difference in cycle period is positive, and reaches a peak of about 
12% midway through the transition. 



Fig. 19.— Normalized difference in cycle-lengths in North and South as function of cycle 
number during change in dipolar to quadrupolar parity. 

This switching is possible because with the circulation pattern chosen, there is much less 
linking of flux between hemispheres at the equator. With dipole symmetry there is strong 
diffusion across the equator, whereas with quadrupole symmetry there is much less since 
both sides have the same sign of held there. Furthermore, dipole symmetry is best preserved 
when the toroidal helds of both signs are strongest where they are being brought into close 
proximity near the equator. This is guaranteed when there is single cell in depth in low 
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and mid latitudes, with equatorward flow near the bottom where the turbulent magnetic 
diffusion is smaller. 


The parity switch ing example we have shown above is very different from that found in 


Jouve fc Brun (120071 1. particularly in latitudes equatorward of 45°. This is not surprising, 
because there is a significant structure in the flow structures - we have two strong primary 
cells at the top and bottom, extending from th e equator to 60° latitu des, associated with two 
weak secondary cells at polar regions, whereas IJouve fc Brun (120071) had four equally strong 
cells (two in top and two in bottom) of equal latitudinal extents. We found a simple drift of 


phase in one hemisphere with respect to the other at all latitudes, while in iJouve &: Rrnn 


(120071) (Figure 10) there is a radical rearranging of poloidal flux in low latitudes, with a 
switch in symmetries by the antisymmetric mode dying out, replaced by the symmetric 
mode without much change in phase locally. This difference suggests there may be multiple 
ways for symmetry switching to occur, which should be explored in the future. 


4. SUMMARY AND CONCLUSIONS 

We have compared flux-transport dynamo model results for five meridional circulation 
patterns that may occur in the solar convection zone, as suggested by solar observations 
and/or hydrodynamic models and full 3D simulations applied to the Sun. We carried out 
simulations for both diffusion and advection dominated regimes. Only the circulation pattern 
is different in each simulation; all other physical processes included are the same. We And 
a wide variety of dynamo behavior, as measured by simulated time-latitude diagrams of 
toroidal and poloidal fields. 

In general, circulation patterns with only one cell in depth and no more than two cells in 
latitude produce the most solar like butterfly diagrams. Two cells in depth leads to antisolar 
butterflies from tachocline toroidal fields, but solar like butterflies at mid-depth where both 
cells have equatorward flow. For this pattern to work for the Sun physical mechanisms 
must exist to inhibit magnetic buoyancy there long enough to allow enough amplification of 
toroidal fields to produce spots, while preventing tachocline toroidal fields from reaching the 
solar surface in any observable form. Four cells in latitude leads to some solar-like magnetic 
patterns, but very fast cycle periods compared to the Sun. Surface Doppler measurements 
also do not support the existence of four cells distributed evenly in latitude, though multiple 
cells confined to polar latitudes can not be ruled out. 

All of the solutions we have found retain dipole or solar-like symmetry about the equator 
within the simulation span of 500 years, except the case of circulation pattern with two cells 
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in latitude and and two in depth when the flow speeds in the upper and lower cells differ 
by less than a certain amount. In that case, from a small difference between hemispheres 
starting probably from the numerical truncation, the solution switches to quadrupole type 
within several magnetic cycles and stays there for ever, even when starting from essentially 
dipole symmetry. This switch is achieved simply by one hemisphere temporarily changing 
its period relative to the other until the relative phase changes by one-half cycle, without 
changing the pattern itself in either hemisphere. A milder version of this effect could be partly 
responsible in the Sun for differences in phase between northern and southern hemispheres 
that do not go so far as to switch the dominant symmetry observed, which is dipolar. 

Despite producing significantly different butterfly diagrams for toroidal and poloidal 
fields, our flux transport dynamo simulations with different meridional circulations have 
many properties in common, as revealed by our parameter survey. In almost all cases for all 
parameters chosen, cycle length monotonically declines with increasing circulation amplitude 
and increasing turbulent magnetic diffusivity, but is nearly independent of poloidal source 
amplitude. Maximum fields generated also decline monotonically with increasing diffusivity, 
but increase with poloidal source amplitude. However, changes in circulation amplitude 
produce non-monotonic changes in peak fields for different circulation patterns. 


Throughout the calculation we have fixed the differential rotation profile given by expres¬ 
sions (11) and (12). A differential rotation profile that more accurately captures high-latitude 
pattern beyond 60° and the near-surface shear layer requires a dditional ter r ns. P erforming 
simulations with such a differential rotation profile as given by ISchou et al.l fll998l) (see also 


Dikpati et ah (1200211 ). we did not And changes in dynamo cycle period or in surface radial 
fields, but a small increase (~ 1.07%) in the tachocline toroidal held amplitude with respect 
to that obtained in the present paper. 

There are at least two important effects related to MHD turbulence that we have not in¬ 
cluded in the model we have used that we need to examine in future studies. Both would add 
nonlinearities to the system that are currently beyon d the scope of formulation in this paper . 


One is the so-called ’turbulent pumping’ mechanism fiTobias et ah II1998I: iKapvla et al 112006 


ul: 


Guerrero fc de Gou veia Dal Pin o 1 120081) and references therein) , and the other is diffusiv¬ 


ity ’quenching’ (IGuerrero. Dikpati fc de Gouveia Dal Pino II2009I) . Apparently the turbulent 
pumping might be able to counter the magnetic buoyancy effect to keep the fields in the 
convection zone, especially for the solutions where the mid convection zone toroidal held 
shows more solar-like butterfly pattern for a meridional circulation with two cells in depth. 
This effect should be explored in the future. 
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